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C^l ■ In this paper, we calculate the cross-section for the photoproduction of 

quark-antiquark pairs in the peripheral collision of ultra-relativistic nuclei, 
by treating the color field of the nuclei within the Color Glass Condensate 
model. We find that this cross-section is sensitive to the saturation scale 
Qj that characterizes the model. In particular, the transverse momentum 
spectrum of the produced pairs could be used to measure the properties 
of the color glass condensate. 
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In recent years, there has been a lot of interest in the description of the dis- 
tribution of partons inside a hadron or nucleus (see |Q, g] for a pedagogical 
introduction). An outstanding issue is the problem of the saturation of the 
parton distributions at very small values of the momentum fraction x [|3[ ^, |3[ . 
Indeed, the BFKL jq, m, pi equation predicts an uninterrupted rise of this distri- 
bution as x becomes smaller and smaller. However, this is not compatible with 
unitarity, but rather an artifact of the linear nature of this equation: in other 
words, the BFKL equation becomes inadequate when the phase-space density 
approaches l/a s ||, [j ^0|. This is an effect of the nonlinearity of QCD, and 
happens even at small coupling. 

Many of the recent investigations are based on a model introduced by McLer- 
ran and Venugopalan O, 021 O, which describes the small x gluons inside a 
fast moving nucleus by a classical color field. This approximation is justified 
by the large occupation number of the soft gluon modes. In this model, the 
classical Yang-Mills equation satisfied by the classical color field is driven by 
the current induced by the hard partons. The distribution of hard color sources 



inside the nucleus is described by a functional density which is taken to be 
Gaussian in the simplest forms of this model. This physical picture describes 
what can be called a "Color Glass Condensate" . The classical version of this 
model depends only on one parameter, called "saturation scale" and denoted 
Q s , defined as the transverse momentum scale below which saturation effects 
start being important. This scale increases with energy and with the size of the 
nucleus jl], ||, [l4J. Therefore, at energies high enough so that Q s 3> A QCD , the 
coupling constant at the saturation scale is a small parameter and a perturba- 
tive approach can be considered. In that case, "perturbative" means that one 
can work at lowest order in a s , but all orders in the large classical color field 
must be kept in principle. 

Many theoretical improvements have dealt with the quantum corrections to 
this model. When going to smaller values of x, one must integrate out the 
modes comprised between the new scale and the former scale |J. Therefore, the 
functional density describing the distribution of hard sources has been found 
to obey a functional renormalization group equation |15[ uM that pilots its 
variations when one moves the separation scale between the hard and soft scale. 
Developments along those lines can be found in 0, g ||, |o[ |l|, || ||. A 



particularly simple and suggestive form of this equation has been found recently 
[ p![ p5[ , as well as approximate solutions p6J. 

On the phenomenological side, this model has been used in order to study 
the fcrmionic degrees of freedom at small x and compute structure functions like 
F 2 [^7| . It has also been used to study the production of gluons in the collision 
of two color sources by several groups of authors. |2^, ^9| calculated the pro- 
duction of gluons at lowest order in the density of soft gluons. Then, Krasnitz 
and Venugopalan |50[ |3l]] attacked the problem of the gluon production to all 
orders, by solving the classical Yang-Mills equation on a 2-dimensional lattice 
for an SU(2) gauge theory. Recently, this problem has been solved analytically 
for the case of pA collisions p2f , for which the gluon density for one source is 
much smaller than that of the second source. However, all these studies deal 
with the central collision of two objects (AA, pA), and the computed number of 
gluons cannot be observed directly. Instead, several layers of complex processes 
are to be taken into account between the stage where gluons are initially pro- 
duced, and the stage where hadrons are observed in a detector (in particular, 
the hydrodynamical evolution, and the freeze-out of the system). 

For this reason, it would be useful to consider a process which is sensitive 
to the distinctive features of the color glass condensate model, and yet can 
be observed in peripheral collisions. The second condition implies that such 
a process involves electromagnetic interactions in some way, because they are 
long-ranged. Since this process should also be sensitive to the gluon density, 
it has to involve quarks, which couple both to photons and gluons. Therefore, 
the simplest such process is the photoproduction of qq pairs, which has already 
been considered as an interesting probe for the gluon distribution [B3l p4| . The 
simplest diagram for this process has been represented in figure [fl. 

In this paper, we calculate this process for the peripheral collision of two large 
ultra-relativistic nuclei. We treat the electromagnetic interaction to lowest order 



Figure 1: Prototype of the diagrams contributing to the photoproduction of a 
qq pair in AA collisions. 



(or, equivalently, to the leading logarithmic approximation), and the interactions 
with the color field of the nucleus to all orders in the classical color field. We find 
that this process can indeed take place at large impact parameters, and is very 
sensitive to the features of the gluon distribution, for quarks whose mass is close 
to the saturation scale. We predict that the transverse momentum spectrum 
da T /dydk± of the produced pairs has a maximum for a momentum close to Q s , 
which could therefore be used as a way to measure the saturation scale. 

The structure of this paper is as follows. In section g, we outline briefly the 
main aspects of the color glass condensate model. In particular, we recall how 
the use of the covariant gauge can simplify the description of the system at the 
classical level. 

In section [| we detail the main steps of the calculation of the cross-section 
for the photoproduction of qq pairs. In particular, we explain how this quantity 
relates to Green's functions evaluated in the classical field, how one can cal- 
culate the retarded quark propagator in this background field, and the impact 
parameter dependence of the average multiplicity. 

Section explains how the transverse momentum spectrum of the produced 
pairs can be used as a probe of the color glass condensate, and section || deals 
with the integrated cross-section. 

Finally, section g is devoted to concluding remarks. 

2 The color glass condensate 

A few years ago, it was suggested that the color field of a strongly boosted 
nucleus can be described as a classical color field |y], [12], [L3| , thanks to the fact 
that the small x gluonic modes have a very large occupation number. 

The classical field A^ describing the soft modes satisfies a classical Yang- 
Mills equation: 

[D li ,F^]=gJ v , (1) 

driven by a source J v generated by the hard modes. Note that this current lives 
in the adjoint representation of SU(N C ): J v = JaT a . For a nucleus moving 
at very high velocity in the positive z direction, this source has the Lorentz 



structured J'' = 5^ + J + . Since the current J^ is covariantly conserved, we have: 

[iW]=0, (3) 

which requiresn 

J + (x + ,x~,xx) = W(x + ,x~,xj_)p(x~ ,x±)W^(x + , x~,x±) , (4) 

with 

p(x~,X±) = J + (x+,x~,Xj_) , 

W(x + ,x~,x±)=Texp<ig dz + A~(z + ,x~ ,x±) \ . (5) 

Therefore, unlike in QED, the color source J + cannot be fixed independently 
of the field A» itself. 

However, it is possible to simplify the problem by using the light-cone gauge 
defined by A + = 0, and by assuming that the physics is invariant under trans- 
lations along the x + axis (d~ = 0) because this is the direction of the nucleus 
trajectory. Under those assumptions, one can find solutions of the Yang-Mills 
equation for which A~ = (and hence W = 1). More explicitly, the remaining 
Yang-Mills equations become 

[Di,F^] -0, 

[D i ,d+A i ]=-gJ + . (6) 

The first of these equations indicates that the transverse components are a pure 
gauge field. Therefore, there is a unitary matrix V such that: 

A l = -V\x^ 1 x^)d l V{x-,x 1 _) . (7) 

By plugging this A 1 into the second equation, one obtains an equation for V in 
terms of the source p. It is however simpler to perform a gauge transformation 
generated by this V: 

A» -»■ 2^ = - Vd^V* + VA^ , (8) 



1 In the following, we make extensive use of the light-cone coordinates. For any 4- vector 
x^ , we define: 

+ x° ± x 3 

x =^T' (2) 

and denote Xj_ the transverse component of its 3- momentum X. With these notations, the 
invariant norm of x M is x 2 = 2x + x~ — X± 2 , and the scalar product of fc M and x 11 is k ■ x = 
k + x~ + k~x + — fcx • X±. The invariant measure d 4 x becomes d^x = dx + dx~d 2 Xj_. 
2 Like for J v , we denote p = p a T a . 
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(9) 



and to the Yang- Mills equation 



-ViA+= 5 p, (10) 

where p = VpV^ is the source in the new gauge. This equation has the simple 
structure of a 2-dimensional Poisson equation with a source p, but remains very 
difficult to solve in terms of the light-cone gauge source p. However, for reasons 
that will become apparent later, it is sufficient to express V in terms of the 
covariant gauge source p: 

V(x~,x ± ) = Tcxp I -ig 2 J dz~ ^paiz' , Xl _)T a \ . (11) 

Having a formal solution of the classical Yang-Mills equation for a given 
distribution p of the hard color sources, the method to use this model is the 
following. One first expresses the observable of interest in terms of the classical 
color field, or equivalently in terms of the covariant gauge source p: 0[p\. Then 
one must perform an ensemble averagen over the set of sources p. The observed 
value of this observable is therefore 

(O) ee J[dp\w[p\0[p\ . (12) 

A simple ansatz for the distribution of color sources w[p\ is to take a Gaussian 
distribution: 

w[p\ = ex P < - dx d x± 2 2 , _, >■ (13) 

In the above formula, p 2 (x~)dx~ has the interpretation of a density of hard 
color sources per unit of transverse area, in the slice between x~ and x~ + dx~ . 
For a very ultra-relativistic nucleus, this function is strongly peaked around 
x~ = 0. At this stage, it is easy to justify the use of either the light-cone gauge 
source p or the covariant gauge source p in order to express observables. Indeed, 
the weight w[p\ is invariant under the gauge transformation of the source, as is 
the measure [dp\. The saturation scale Q 2 S is related to the integral of p 2 (x~) 
over x~ (see Eq. (70) in the appendix W for the definition). Its value is expected 



to be of the order of lGeV at RHIC and 2 - 3GeV at LHC. 



3 In a given collision, the hard color source appears to be frozen due to time dilation. 
However, this frozen configuration is going to change from event to event. 



With such a weight, averages of products of sources p can all be expressed 
in terms of the average of the product of two of them 



(p a (x ,xj_)p b (y ,y±)) = 5 ab p 2 (x )S(x -y )S(x±-y±) 



(14) 



thanks to Wick's theorem. More precisely, the average of the product of an even 
number of p's is obtained by performing the sum over all the possible pairings 
of two jcTs, and the average of the product of an odd number of /J's is vanishing. 



3 Calculation of da T /dy 



3.1 Relation to Green's functions in the classical field 

Let us start by giving the connection between the average number of pairs pro- 
duced per collision in terms of a fermionic correlator in the classical background 
field. For a given configuration of the hard color sources, one can use the for- 
malism developed in |3q| to express this quantity in terms of the retarded quark 
propagator. Using the same notations as in |35||, we have 



i\p\ = 



d 3 q 



d 3 p 



(27r) 3 2w g J (2tt) 3 2w p 



u(q)T<p(q,-p)v(p) 



(15) 



for the average number of pairs produced in a collision at impact parameter b. 
In this formula, Tjf denotes the interacting part of the retarded propagator of 
a quark in the classical field generated by the covariant gauge color source p and 
the electromagnetic current of the other nucleus. Note also that this quantity 
can be made differential with respect to the quark or the antiquark momentum 
by not integrating over the corresponding momentum. This will be useful in 
order to calculate the multiplicity per unit of rapidity. 

Then, averaging over the configurations p of the hard color source, one ob- 
tains n: 

' k/H „<[,->] "[/'] • (16) 

Note that n depends on the impact parameter b of the collision. By inte- 
grating over impact parameter, one obtains the total cross-section. 



3.2 Retarded quark propagator 

In order to calculate n, we need the retarded propagator of a quark moving in 
a classical field made of the electromagnetic field of one nucleus, and the color 
field of the other nucleus. In |3q] , we have shown that the retarded propagator 
becomes very simple in the ultra-relativistic limit. This result was derived for 
QED, but since the argument was based only on causality, it can be extended 
immediately to QCD. Therefore, we can write formally (convolution products 
are not written explicitly) for the scattering matrix Tr: 



1 R - 1 R 



T 



QED 



1 R ^R 1 R 



QED n QCD 

R. ^Rl-R 



(17) 



where T R and T R are the retarded scattering matrices of the quark in the 
electromagnetic field of the first nucleus and in the color field of the second 
nucleus respectively. G R is the free retarded propagator of a quark, which has 
the following expression in Fourier space: 

g °r(p) = i 2 ^ + 2 T 9 o • ( 18 ) 

Assuming that the nucleus acting via its electromagnetic field is moving along 
the negative z axis, we have in the ultra-relativistic limitn 

T R ED {q,p) = 2ttS{ P + - q+)ij- f d 2 xAe- leAt - x ^ - l\ e ^±-P±)^± , (19) 

where A(x±) is the 2-dimensional Coulomb potential of the first nucleus. Note 
that this potential can also be written as 

A(au.) = I dz-^ jT p{z-,z x _), (20) 

J-oo V _l 

where p(z~ , z±) is the electric charge density inside the nucleus. 

Since the second nucleus is also ultra-relativistic, we can assume that the 
quark does not resolve the longitudinal structure of the color source. The only 
difference with QED in that respect is the ^"-ordering. If one uses the gauge de- 
fined by equations (|9|) , then the problem of the propagation of a quark through a 
color glass condensate is formally similar to the QED case, and we can generalize 
equation (|l^) without difficulty: 

T R CD (q,p) = 2nS(p~-q-)J + Jd 2 x ± [u(x ± )-l\e^-i>^ , (21) 



with 



U(x±) =Tcxpj -ig 2 / dz ^-paiz ,x ± )t a i 



(22) 



This formula also takes into account the fact that this nucleus is moving in 
the positive z direction. The inverse Laplacian acting on the covariant gauge 
color source, integrated along the axis of the collision, is the analogue of the 
2-dimensional Coulomb potential in QED (see Eq. (EQ)). Note also that the 
covariant gauge color source is now contracted with the generators of the fun- 
damental representation. More details on the eikonal approximation in a color 
glass condensate can be found in p6[ . 

For pair production, only the last two terms in Eq. ( |17|) are relevant, because 
one cannot produce a qq pair on-shcll by interacting with only one nucleus. After 
doing the integrals over the + and — components of the momentum transfer 

4 We denote v± = (1, 0, 0, ±l)/v / 2 the 4-velocities of the two nuclei. 



(made trivial by the delta functions) , the relevant part of the retarded scattering 
matrix isR: 

T^{q, -p) = -ij j^fF R {k±)F R (p ± + q±- k±) 

l2g-p+ + (q_L - fcjj 2 + m 2 2q+p- + (p±_ - fcj_) 2 + m 2 / 
where we denote 



F R CD (k ± ) = Jd 2 x ± [U(x ± ) - 1] e lk 



OED 



F R {k±)= d 2 x ± e _ * eAla! - L ' - 1 e**- L - iB - L . (24) 

In the following, it will prove convenient to factor out the impact parameter 
dependence contained in F R (we chose the origin of the transverse coordinates 
in such a way that the nucleus acting via its color field has its center at x± = 0, 
so that all the b dependence goes into F R ), by writing 

i^ D (fcx) = e^Vr D (fcx). (25) 

Note that, at the Born level, we have 



n QED 



AirZa 

Born fc 2 



fn ( fc ^)„ =-i3-> ( 26 ) 



where a — e 2 /47r is the electromagnetic structure constant. The infrared 
sensitive denominator l/fej_ produces a large logarithm^ in the cross-section. 
Coulomb corrections beyond this Born term do not produce such a large loga- 
rithm J37], pq, p9], and will be discarded in this paper. 

3.3 Expression of n 

Using the formula of Eq. ([l6|), we can write 

d 3 p d 3 q f d 2 k± f d 2 l± J „qcd *qcd 



f d A p d 4 q f d z k ± f d 2 l ± / ecu > pt QCD /, n 

J wtm> Wn^ JWrJ w? TTc \ Fr {kl - )FR (z±) 

QED QED 

xF R (p±+q±-k ± )F R (p± + q_L-l±) 

xTr (M(p,q\k ± )M*(p,q\l ± )), (27) 



5 We have changed p — ► — p as required in formula Eq. (P-3). In addition, we have also 
changed k± —> q± — k± in the first term, and k± —* Px + &x i n the second term for 
notational convenience. 

6 By "large logarithm", we mean a logarithm containing the center of mass energy of the 
collision s, or the corresponding Lorentz factor 7. 



where Tr c is a trace over color indices, Tr is the Dirac trace, and where we 
denote 



M(p,q\k ± )=u{q){ 



y+(«/-#4 


m,y_ 




2q~p+ + (gjL - 
y_(jt- 


fc^) 2 + 771 2 



2< Z +p- + (p ± -fc ± ) 2 + m 2 H P) - (28) 
Using the results of appendix [A] (see Eqs. ( |74| ) and (|S0"|)), we can rewrite this as 
<i 3 p d 3 q /" d 2 kj 



N 
' ' (27r) 3 2w p (27r) 3 2w g J (2tt) 2 J (2tt) 

x f d 2 x ± d 2 y ± e lk ^ x ^e- a± - v± P(x ± )P( yi _) [l + e -^(^-^) _ 2e 

xFr (p_l + q_l - fej.)-FR (p_L + g_L - Jj_) 
xTr(M(p, Q |fe ± )M*(p )Q lZj.)) , (29) 

where the function V(x±) describes the transverse profile of the nucleus (see 
the appendix |A| for more details on the effects of the nucleus finite size) . For the 
terms in the square bracket which are independent of x± and y±, the integral 
over those variables is trivial and just leads to the product V(k±)V(l±) of the 
Fourier transforms of the functions V . For a large nucleus, the Fourier transform 
V(k±) is very peaked around k± = 0, with a typical width of 1/R, where R is 
the radius of the nucleus. However, it is straightforward to check that 

M(p, q\0) = M(p, q\p^ +qx)=0. (30) 

Moreover, the only scale controlling the amplitude M(p, <7|fcj_) is the mass m of 
the produced quark, which means that it varies significantly only if k± changes 
by an amount comparable to m. Therefore, for heavy quarks such thatn m 3> 
1/.R, the function M(p,q\k±) remains very close to zero over all the range 
over which V(k±) differs from zero. As a consequence, we can neglect the 
contribution of the terms independent of xj_ and y±, and focus on the term 
proportional to exp(— B2(x± — y±_)). 

Using now the approximation of Eq. ( J76| ) for the function^] B2(x± —y±), we 
see that this function becomes much larger than 1 as soon as \x± — y±\ 3> Qj 1 - 
Since we have also R ;» A~* D ;» Qj 1 , it is reasonable to approximate the 
product V{x±)V(y±) by a single factor V(x±). The integral over x± and y± 
then becomes trivial and gives 

U Nc J ( 27 r) 3 2u, p (2^) 3 2u, g J (2^) 2 J (2^ V{k± ' ±)C(fe±) 

QED QED 

x/r {P± + q±-k±)f R (p± + q±-l±) 
xTk(M(p,q|fc ± )M*(p > q|I ± )) ) (31) 



7 Note that for a large nucleus of radius R = 6fm, one has 1/2R ~ 17MeV. 
8 For this approximation to be valid up to distances \X± — y±\ ~ Q7 , one needs Q s 2> 
^■qco- 



where we have defined 



C{k x )= J d 2 x ± e ik±x± e- B ^ x ^ . (32) 



Using the Eqs. (||), @), (JtTJ) and @ of the appendix 0, one can see that 
C(k±) is nothing but 

C(k±)= I d 2 x ± e ikx - xx (U(0)U*(x±)) . (33) 

In other words, by assuming that the nucleus is large (R 3> A~* D ), we have been 
able to factorize the dependence on the nucleus size and recover the correlator 
(U(0)U'(x±)) associated to an hypothetical infinite nucleus. 

3.4 Impact parameter dependence 

From Eq. (pl|), one can determine the main features of the impact parameter 
dependence of the multiplicity n. It is useful to recall that M(p, q\k±) vanishes 
when p± + qj_ = k± . Moreover, since the typical transverse momentum p± + 
q± — k± has to be very small for coherent photonsg, compared to the other 
scales in M(p, q\k±) which are controlled by the quark mass, on can use a 
Taylor expansion of this amplitude around the zero: 



M(p, q\k±) « (pjl +q±- k±) ■ L{p, q) , (34) 

L(p,q) = u(q) ' 



with 



2(p+ + q + ) V T V- 

In addition, the small transverse momentum of the photons allows one to write 
C(k±) w C(p± + q±). Therefore, one can see that the impact parameter de- 
pendence of n is completely controlled by the integral 

[) -J (2tt)2 (2^)2 6 ' (mi - n±) ml n\ 

(36) 
where we have used the new variables m± = p± + q± — k± and n± = p± + 
q± — lj_. If the impact parameter is very large (6 > i?), then the momentum 
difference \m± — n±\ can at most be of order 1/6 <C 1/R because otherwise 
the exponential would average to zero. On such a small range, the function 

9 In order for a photon to be coupled coherently to the total electric charge of the nucleus, 
its transverse momentum should not be larger than 1/2R. For R = 6fm, 1/2R ~ 17MeV. 
Note that the lower bound on the transverse momentum of the photon is 771/7, where 771 is 
the mass of the produced particle. 



10 



V(m± —n±) can be replaced by "P(O) = J d 2 x±V(x±) = irR 2 , and the integrals 
over mi and n± are trivial to perform. One obtains: 

Tlh\ -T>2 bL (P,<l) b-L*(p,q) 

1{b) = * R ~^b 2 2^ 2 ~ ■ (3?) 

Therefore, the multiplicity decreases like 1/b 2 at large impact parameter. If 
taken seriously up to b = +00, this would lead to a logarithmic divergence in 
the cross-section after integration over fo, as is well known in problems involving 
long ranged electromagnetic interactions. In reality the field does not keep its 
Weizsacker- Williams form up to arbitrarily large transverse distances. Indeed, 
at the distance b from the source, the field has a typical longitudinal extension 
of 6/7, where 7 is the Lorentz factor of the electromagnetic source. This width 
should be compared with the typical Compton wavelength of the particle being 
produced, i.e. 1/m in the present case, which leads to an upper cutoff b = 7/771 
in the integral over impact parameter. This upper bound on b is equivalent to 
the lower bound m/7 on the transverse momentum of the photon. On the lower 
side, one must cutoff the integral at b = 2R in order to consider only peripheral 
collisions. Although these cutoffs are introduced by hand, this ansatz gives the 
leading logarithm correctly. This leading log approximation is valid if there is a 
large range between the scales 2R and 7/777. At smaller impact parameters, i.e. 
when b is close to 2R, the 1/b 2 scaling law is modified by edge effects which are 
difficult to estimate analytically. 

3.5 Evaluation of da T /dy 

In order to obtain the leading logarithm of the cross-section for peripheral col- 
lisions (the ones that are phenomenologically interesting for photoproduction), 
we can use the formula Eq. (g7|) for all the range 2R < b < 7/777. Noting that 

J d 2 b VVf(b 2 ) = S -^Jd 2 b b 2 f(b 2 ) , (38) 

we can write 

• b L b-L* L L* f d 2 b 

Therefore, we have the following expression for the totaf^ cross-section for qq 
photoproduction in peripheral collisions: 

7/m 7/m 

d 2 b 



g t = J d 2 bn = irR 2 2N c (Za) 2 f 

2R 2R 



10 By "total", we mean the cross-section obtained by counting all the produced pairs. This 
quantity in general differs from the cross-section obtained by counting only events where 
exactly one pair is produced. However, when the average multiplicity is small, the two are 
very close, and this distinction becomes pointless. 
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This is the formula we are going to use in the following. In Eq. (pQ), the factor 
which is the less well known analytically is the correlator C(p±+q±). Therefore, 
in the evaluation of integrals, we are going to keep the variable p± + q± (i.e. 
the total transverse momentum of the pair, or approximatively the momentum 
transfer between the color field and the qq pair) unintegrated until the end. 
At this stage, we will have the choice to use some analytic approximation of 
C(p_l + Q±) or do the integral numerically. This approach also leaves open the 
possibility to consider other saturation models. 

The first step is to perform the Dirac algebra, which leads to: 



Tr (\L(p,q) 



(p q 

p~q~(p± ^ 



( P + + q+)y 



q±) 2 -2 



(q+q p x +p + p q±) 



(p-* 



7+^2 



(41) 



One can check that the expression of Eq. (|d]) vanishes for p± + q±_ =0, i.e. 
when the momentum transfer from the color glass condensate is zero. By using 
the fact that the final state particles are on-shell, we can perform for free the 
integrals over q + and p + . Then, the integral over p~ is very easy because the 
dependence on this variable is rational. Noticing finally that dq~ /q~ = 2dy, 
where 2y = ln(g + /g _ ) is the rapidity of the produced quark, it is very easy to 
obtain the differential cross-section with respect to the rapidity 11 of the quark 
(or equivalently, of the antiquark) : 



da T 
dy 



ttR 



7/m 

, 8N c (Za) 2 r fb 
V 2 J ~~& 

%R 

„ i (P± + <lx? 



d 2 p ± d 2 qj_ 
17^S C fa 



H"l 



2tt) 
Pi 



(27T) 

91 



■qx) 



(42) 



with loz 



Pi 



m 2 and u> 2 



Let us take as one of the integration 



variables the total transverse momentum of the pair k± = p± + q± , and keep 
this variable unintegrated until the end. Then, it is easy to perform the integral 
over the vector q± . An intermediate step is the integral over the angle 9 between 
q± and k±: 




fe 



M-X 



2 
Tx 



fei 



2k 2 - 



6(g 2 + m 2 )^(m 2 + q 2 +k 2 ) 2 - 4q 2 k 2 



11 The rapidity distribution is found to be flat here, because we do not impose any upper 
bound on the energy of the quark. However, by arguing that the produced particles cannot 
be more energetic than the incoming nuclei, one of course finds that the rapidity distribution 
cannot extend beyond the rapidity of the projectiles, i.e. iln('y). Integration over rapidity 
would therefore bring another factor 111(7) i n the total cross-section. 
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m 2 (,i + fci) + (,i-fci) 2 to 2 (43) 



Integrating then over the modulus q±, we obtain our final expression for da T /dy: 



da T 
dy 




o2 ^v c (W 


7 /m 

/• d 2 b 
7 "b 2 " 








27? 








+ 00 


I 


f 




X 


/ dk ± k ± C(k_ 



J l 


1 + 



4(fc 2 -m 2 ) , jfcj. I ,^ 

7 =arcth } .(44) 



k± \Jk\ + Am 2 y/k'j_ + 4m 2 

At this stage, it may be useful to summarize all the assumptions concerning 
the various parameters that entered in the derivation of this result. First, we 
have assumed that the nucleus is large 

R»Kco (45) 

in the average over the color sources (appendix [A|) and before Eq. (|l|). This a 
very good approximation for a nucleus like gold or lead. The next assumption 
is that 

Qs » A QCD , (46) 

which has been used in the approximation that lead to Eq. (|3l|). We take 
A QCO w 0.2GeV, while Q s is expected to be of the order of lGeV at RHIC 
and 2 — 3GeV at LHC. Therefore, this assumption is only marginally satisfied 
at RHIC, and probably much better at LHC. We have also assumed that the 
quark mass is large enough 

m > R- 1 (47) 

before Eqs. (M) and (f40|). In addition, for the leading log approximation in the 
integral over impact parameter to be justified, one needs a large enough center 
of mass energy 

7 > mR . (48) 

4 k± spectrum 

In Eq. (144), da T / dy is given by a 1-dimensional integral involving the correlator 
C(k±) defined in Eqs. (|32| ) and (p3|). If all one needs is a spectrum as a function 
of the total momentum of the pair, then it is not necessary to integrate this 
expression further, and the k± spectrum is directly given by 



Tri? 2 —°v—/ l n (JLJ) 
\2mRJ 



da T __ Da 2N e {Za)_ 

m l2™ 

TO 2 ) 
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arcth - } . (49) 

y/kl + 4to 2 J V ' 



We observe that this spectrum has the interesting property of being directly 
proportional to the correlator C(k±). According to the discussion of appendix 
pi this correlator is sensitive to saturation effects as soon as k±<Q s . The 
possibility to observe saturation effects by measuring the k± spectrum of qq 
pairs produced via photoproduction is illustrated in the figure @, which contains 
plots of the k± spectrum for the charm and bottom quarks. One can see a 
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Figure 2: Plot of da T /dydk± as a function of k± /A QCD , in arbitrary units (only 
the k± dependent factors of Eq. (49) have been plotted). The value of Q s /A QCn 
is set to 10. The dotted curves are based on a numerical evaluation of the full 
C(k±) as predicted in the color glass condensate model. The continuous curves 
have been obtained with the asymptotic form C(k±) = 2Q 2 s /k\ : which is the 
lowest order in Q 2 S . Open dots and solid curve: m/A QCD = 23 (bottom), filled 



dots and dashed curve: to/A 



QCD 



8 (charm), with A 



QCD 



0.2GeV. 



qualitatively important difference between the prediction based on the color 
glass condensate (dots) and the result one obtains by assuming no saturation 
(solid curves), i.e. by using only the first order term in Q 2 S for C(k±): 



C(k ± )^2^-, 



(50) 



derived in appendix [b|. The color glass condensate model predicts a maximum 
in the k± spectrum, for a value of k± closely related to the saturation scale Q s . 
It is in fact easy to understand why, for a large mass to, the location of the 
maximum is completely controlled by Q s . Indeed, if the mass is very large, the 
function that appears inside the curly brackets in Eq. (149) can be approximated 
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by: 



4(k 2 ± - to 2 ) 



: arcth ■ 



&-L\Ai +4to 2 ^fc 2 ^ + 4to 2 



7k\ 



k±<^m 



D TO^ 



fei 



(51) 



This indicates that the mass dependence becomes a trivial factor m~ 2 which 
does not affect the location of the maximum. The location of the maximum is 
entirely determined by the function kj_C(k±), and therefore depends only on 
Q s . A semi-analytic estimate of the value of k± at the maximum leads to 



fc max 



Ws-G£ 



(52) 



For intermediate masses, this is not exactly true, but we have checked numer- 
ically that the location of the maximum varies very little. For instance, for 
Q S /A QCD = 10, the maximum varies from kj_/A QCD = 9 to k±/ A QCD = 12 if 
one varies the mass between m/A QCD — 5 and m/A QCD = 100. Observing ex- 
perimentally this maximum could therefore be an evidence for saturation, and 
provide a measurement of the saturation scale. 



5 Integrated cross-section 

In this section, we discuss various properties of the cross-section integrated over 
the transverse momentum of the pair. In particular, we provide analytical results 
for the limit of large mass (to 3> Q s ) and the limit of small mass (m < Q s ). 
To that effect, we will need also the large k± expression of the function that 
appears inside the curly brackets in Eq. (|44|): 



4(fc 2 



rcth ■ 



kjL\fk\ +4to 2 y/k\_ + 4to 2 



41n(Jbj./m) 



(53) 



k±^>m 



5.1 Large mass limit 

In the limit of large quark mass, the typical k± in the integral of Eq. ( B ) is 
of order m or smaller. In addition, since m ^> Q s , then so will be k±. For 
such large values of k±, we can use Eq. ( J50| ) for C(k±). Note also that C(k±) is 
strongly suppressed compared to this asymptotic value if k± < A QCD (see figure 
in appendix O), which provides a natural lower bound for the integral over 
fej.. 

Making use of Eq. ( pl| ) and performing the integral over k± , we obtain the 
following asymptotic expression for da T /dy when m^> Q s : 



da T 
dy 



= nR- 



m>Q s 



UN c (Za) 2 Ql 

97T 3 TO 2 



hi 



A, 



In 



2mR 



(54) 
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Note that this formula cannot be used for very large masses, such that the 
argument j/2mR of the logarithm becomes close to 1 (or smaller) . Indeed, the 
logarithmic approximation used in this paper is no longer valid. Moreover, the 
kinematical lower bound 771/7 f° r the transverse momentum of the photon being 
larger than the upper limit 1/2R of the Weizsacker- Williams photon spectrum, 
the photons that can produce such a massive pair of particles has to come from 
the tail of the photon spectrum, which is exponentially suppressed. Details on 
this can be found in [[fO) where is studied the photoproduction of top quark 
pairs in peripheral heavy ion collisions. 

5.2 Small mass limit 

Similarly, one can deal analytically with the case of a very small mass (to <Qs)- 
For this case, it turns out to be simpler to revert the integral in Eq. ( [J4|) to 
transverse coordinate spacepl 



-'*^"( 5 2 5 W*-('--^ , )( J ir) 



dy 



(55) 
where the last factor is nothing but the inverse Fourier transform of Eq. (p3J), in 
the range x± <C 1/m. For x± 3> 1/w, on can check that this Fourier transform 
is strongly suppressed, hence the phcnomenological upper bound in the integral 
over x±. Because of the behavior of the function B2(x±), the combination 
1 — exp(— B2(x±)) behaves mostly like 6(x± — 1/Q S ) where 9 is the Heaviside 
step function. Therefore, we find for the small mass limit: 



da T 
dy 



m<<Q 3tt 2 \mj \2mR 



H 



=^^™V%W^). w 



One can note that the integrated cross-section could also have been obtained 
from the expression of the structure function F 2 (x,Q 2 ) derived in J27|]. Taking 
the large and small mass limit in F2, and convoluting with the Weizsacker- 
Williams spectrum, one can recover the Eqs. ( p4J ) and (pq). 

5.3 Numerical results 

One can also evaluate numerically the integral over transverse momentum in 
Eq. ( B ) . The result of this calculation is plotted in figure g. One can see that 
the asymptotic formula obtained in Eq. ( pi) ) is reasonably accurate for to > Q s , 
up to an additive constant correction to the \n(m/A QCD ) which could not be 
determined analytically and has been fitted. 

In addition, one can see that for large nuclei (R = 6fm) at LHC energies 
(7 = 3000), the value of the cross-section is about 800mb for charm quarks, and 

12 Because the expression between the curly brackets vanishes at fcj^ = 0, it is possible to 
replace e~ B2 by e~ B2 — 1 in the definition of the correlator C(fcj_). We have made use of this 
freedom here. 
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Figure 3: Plot oida T /dy as a function of m/A QCD . In this plot, Q S /A QCD = 10. 
We use R = 6fm and 7 = 3000. The open dots are a numerical evaluation of the 
integral in Eq. (Q) , while the solid curve is the asymptotic formula of Eq. (J54|) , 
in which we have replaced ln(m/A QCD ) by \n(m/A QCO ) + 0.92 (the additive 
constant has been adjusted so that it fits at large values of to). We have circled 
the points corresponding to the mass of the charm and bottom quarks, assuming 



A, 



0.2GeV. 



lOOmb for bottom quarks. This is before any experimental cut is applied on the 
range of k± . 

6 Conclusion 

In this paper, we have evaluated the cross-section for the photoproduction of 
qq pairs in the peripheral collision of two ultra-relativistic heavy nuclei. Both 
nuclei are treated as classical sources, respectively of electromagnetic field and 
color field. The classical color field is described using the color glass condensate 
model. 

It appears that this process is very sensitive to the properties of the gluon 
distribution in the region where saturation might play an important role. In 
particular, we find that the saturation of the gluon density modifies significantly 
the shape of the k± spectrum of the produced pairs. We predict a maximum in 
da T /dydk± when k± is comparable to the saturation scale Q s . Therefore, one 
could possibly use this process as a test of the color glass condensate model, 
and as a way to measure the saturation scale. 
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It is also interesting to study the case of the diffractive production of qq pairs, 
where the exchange between the color source and the pair is color singlet and 
carries no net transverse momentum. One would then have events where the 
quark and antiquark have opposite transverse momenta, with no other object 
in the final state. This calculation is a work in progress. 
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A Calculation of (U(x ± )W( y± )) 
A.l (U( x± )} 

In this appendix, we calculate the correlator {U(x±)W(y)) with 

^+00 



U(x±) = Tcxp I -ig 2 / dz ^p a (z ,a:±)i a ^, 



(57) 



where t a is a color matrix in the fundamental representation of SU{N C ). More 
explicitly: 



— T Pa(z ,x±)= d 2 z±G a 



(x± - z±)p a (z ,z±) , (58) 



where Go is the propagator associated to the 2-dimensional Laplacian. This 
propagator satisfies: 

d 2 
—^-G (x± - z_l) = S(x± - z±) . (59) 

More explicitly, we have 

G (x ± -z ± ) = -Jj^ ^ . (60) 

Let us first consider the average of U(x±). It is in fact useful to calculate 
the slightly more general quantity (U(a~, b~\x±)), where we define 

U (a-, b- \x ± )= Tcxp I -ig 2 f dz' ^p a { Z - , Xl _)t a \ . (61) 

By expanding the time-ordered exponential, we have 

(U(a-,b-\x ± )) = Y, iZ! T^ / IT [d 2 Zi±G Q (x ± - z t± )] 

n=0 ' •* i=l 
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dz-. 



dz 



dz n (p ai {z 1 , ZlJ.) • • • p an (z n , z n ±)) t ai ■■■ t a " . 

(62) 



The bracket under the integral can be reduced to averages of products of two pi's 
thanks to Wick's theorem. To make the discussion more visual, we represent 
the Wilson line integral in U by a bold straight line, the free 2-dimensional 
propagators Go by wavy lines, and the elementary correlator (p a pb) by a cross 
inserted between two propagators: 



sT®-^ 



> £ > 



= G {X± - Z 1± )G (X± - Z 2 ±) (paA z l , z l±)pa 2 {z 2 , «2±)) 



Zl 



<-2 



(63) 

Because of the S(z 1 — z 2 ) contained in such a contraction, we will call it a 
"tadpole" in the following. The only possibility for multiple pairings on a line 
ordered in z~ is to have adjacent tadpoles, like in the following diagram: 



4*-^ 



.^8h> J; 



xS^T& 



Z7 



Z2 Z 3 



Z 4 Z 5 



2+ 



Z~6 



Diagrams with nested or overlapping loops are vanishing because they have no 
support on the z~ axis. Noticing also that 



dz 2 S(z 1 -z 2 ) = 



1 



(64) 



we finally obtain 
(C/(a _ , 6~|:r^)) = exp ■ 



-(*X) 



dz p (z ) 



d 2 z ± Gl{x 1 _ -z±) 



(65) 

We can recognize in this expression the color density per unit of transverse area 
obtained by integrating p 2 along the z~ axis between a~ and b~ . 

A.2 (U(x ± )W(y ± )) 

Let us now turn to the case of the average of the product U (x ±)W (y ±) . One can 
have pairings that connect the Wilson lines corresponding to the two IPs, but the 
z~ ordering prevents the corresponding links to be crossed. We therefore have 
only "ladder" topologies, or "tadpole" topologies similar to those encountered 
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above. In summary, wc need to rcsum all the topologies like: 



•* 4 f C', 



Diagrams that do not fall in this class vanish because of the ordering. In order 
to evaluate the correlator (jjW), we proceed as above, noticing that for a link 
connecting the two Wilson lines, the integral over z~ does not produce a factor 
1/2 because there is no relative order between the z~'s on the two lines. Then, 
we can write the correlator as a sum over the number n of rungs in the ladder: 



(U(a-, b-\x ± )U\a~, b~\y±)) = 



71=0 



dz l 



dz~ 



x(U(a ,z 1 \x±))(U(z 1 ,z 2 \x±))---(U(z n ,b \x±)) 
x (U(a~ , Zi\y±)) (Ufa , z 2 \y±)) • • • (U{z~ ,b~\y±)) 



g\t a t a ) I d z z ± G ( Xl _ - z ± )G (y± - z±) 



fJ- Oi ) • ' ' M {z n ) 



(66) 



where we have already resummed all the tadpole insertions on the Wilson lines 
between the rungs of the ladder. One can notice that all the color matrices come 
in pairs t a t a , which are proportional to the unit matrix. Due to their exponential 
structure, the tadpole insertions can be combined into a single exponential which 
depends only on the endpoints a - and b~ , and can therefore be factored out of 
the integral. The sum over n then trivially exponentiates, so that: 

(U(a-,b-\x ± )U\a-,b-\y ± )) = (U(a~, b~\x±)) (U "(a" , b~ \y ± )> 



xexp{g 4 (t a t a ) 



dz /i 2 (z ) 



/ d 2 z ± G a (x 1 _ - z±)G (y± - z A 



(67) 



Combining all the factors, taking a 
(U(x ± )U\y ± )) 



-oo, b 



+oo, we can write: 



--exp{-^(t a t a ) 



+ OC 

fdz-n 2 (z-) fd 2 z ± [Go(xl-zl) - G (y±-z ± )Y 



(68) 
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One usually introduces here the saturation scale Q s , defined byrj 



+ OG 

)l = y(*X) fdz-^(z-) . (70) 



A. 3 Finite size effects 

The above calculation has been performed under the assumption that the color 
source has an infinite transverse extension, which leads to translation invariance 
in the transverse plane for the correlator: 

(U(x±)Ui{y±)) = (U(0)U*(y± - x±)) . (71) 

In order to see how the finite size of the nucleus enters in this correlator, we 
have to go back and replace /j, 2 (z~) by /j, 2 (z~)V(z±), where the function V(zj_) 
describes the transverse profile of the nucleus (if one neglects edge effects, this 
function is 1 inside the nucleus and outside). Because of our ansatz to factor- 
izePjthe profile from the density [i 2 (z~), the transverse profile does not change 
the calculation of the correlator. Simply, the argument of the exponential of 
Eq. (B8h becomes (up to a sign): 



A(x±,y±) =Q 2 S f d 2 z ± P(z ± ) [g (xj_ - z±) - G { yi _ - z±j\ * . (72) 

We see that the finite size of the nucleus reduces the support of this integral 
to the transverse section of the nucleus. Inspection of this integral indicates 
that there is a residual logarithmic infrared divergence (a stronger, quadratic, 
divergence is cancelled thanks to the difference of the two propagators) . This 
divergence is in fact well known and it has been argued that it is naturally 
screened by confinement (a feature which is not present in the color glass con- 
densate model) at a scale comparable to the nucleon size (see 41 for a more 
elaborate discussion of this issue), which is of the order of A~J, D . Practically, 
one can achieve this by cutting off the propagator Gq(x± — z±) in order to make 
it vanish for distances larger than A"* . For instance, one can regularize the 

13 We have chosen here a definition of Q 2 , which is particularly convenient when coupling 
the color field to fermions, because it absorbs the Casimir in the fundamental representation. 
Other conventions have been used in the literature. For instance, Ja defines 

Q 2 = 2irN c a 2 / dz~^(z~). (69) 

— oo 

14 This is of course not exact because for a spherical nucleus the transverse profile depends 
on z~ (or, equivalently, Q s is smaller at the edges of the nucleus than at the center), but this 
ansatz captures the main features due to the nucleus finite size. Indeed, in this particular 
problem, the photon that hits the nucleus has a very small transverse momentum and cannot 
resolve the transverse size of the nucleus. This photon therefore sees a density of color charges 
averaged over the transverse section of the nucleus, or an averaged saturation scale Q s . This 
averaged Q s is dominated by the rather large value at the center of the nucleus. 



21 



Fourier transform in Eq. ( |60[ ) by replacing kj_ by fcj_ + A 2 D , which leads to a 
propagator proportional to Kq(A qcd (x± — z±)). All the predictions made by 
using this prescription will depend on A~\ in a universal way (i.e. indepen- 
dently of the details of the implementation of the cutoff) as long as one looks 
only at momenta large compared to A QCD . 

Under the assumption that the nucleus is very large (i.e. its radius R satisfies 
R ~S> A~* D ), one can split the function A(x±, y±) according to whether x± and 
y± are inside or outside of the nucleus: 

A(x±,y±)= (1-V(x ± ))(l-V(y±))x0 
+P(x x )(l-V(yx))B 1 (x x ) 

+ {l-P{x ± ))P{y ± )B 1 (y ± ) 
+V(x ± )T(y±)B2(x ± - y ± ) , (73) 

where we denote 

Bi(asj.) = Ql I 'd 2 z ± G 2 Q ( Xl _ - z ±) , 

B 2 (x ± ) =Q 2 S J d 2 z ± [G (z ± ) - G Q (z ± - x ± )} 2 . (74) 

In Eq. (J73), the neglected terms are suppressed by powers of (RA QOD )~ 1 , Using 
now the Fourier transform representation of Go, it is trivial to arrive at 

+00 
£ 2 (*J-) = — I %{1~Hpx^)). (75) 

7T J P A 



Cutting off the logarithmic divergence by A QCD , we find 

2M*x)«^W— J— ) . (76) 

4tt \x ± A qcd J 

an approximation valid for transverse separations x± much smaller than A"*^ . 
For B\, we find to lowest order in A QCD /Q S : 

B ^ ~ TT~ » l ■ (?7) 

QCD 

Using those results, we can rewrite the two correlators considered in this 
appendix as: 

(U{x ± )) = (U\x ± )) = 1 - V(x ± ) + T(x ± )e- B i , (78) 

and 

(U(x ± )uHy ± )) = (1 - V(x ± ))(l - 7>(y ± )) 
+T>(x ± )(l-V(y ± ))e- B i 
+(l-V(x ± ))V(y x )e- Bi 
+-p(x±)V(y±)e- B * {x± - y±) . (79) 
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From this, it is trivial to obtain the correlator needed in section 

((U(x ± ) - l)(U*(y±) - 1)) = P(x ± )V(y ± ) [l + e r B ^^^ - 2e~ l 

J (80) 

Note that in the derivation of n, the factor U(x±) — 1 comes from the retarded 
amplitude Tr while the factor W(y±) — 1 comes from its complex conjugate. 
The coordinates x± and y± can be seen as the transverse coordinate of either 
the quark or the antiquark. Therefore, the proportionality of this correlator to 
the factor V(x±)V(y±) simply means that in order to produce a qq pair, either 
the quark or the antiquark has to intercept the nucleus that acts via its color 
field (and that, both in the amplitude and in its complex conjugate). 

B Asymptotic expansion of C(k±) 

The pair multiplicity depends crucially on the correlator defined by 

C(k±) = I d 2 x ± e lx± - k± e~ B2(x±) = I d 2 x ± e lx ^- k ^(U(0)U^{x ± )) . (81) 

This object depends only on the modulus of fcj^ thanks to the rotation invariance 
of the function B^ ■ 

It is rather easy to obtain the first two terms of the asymptotic expansion 
of C(k±) for k± 3> Q s - Indeed, for such a large fcj_, the typical x± contributing 
in the integral is of order l/fcj_, and Q s x±_ <C 1. We are therefore in the region 
where the function B2(xj_) is know analytically through formula (|7q) . The 
asymptotic expansion of C(k±) is then obtained by expanding the exponential 
of — B2, which gives 

One can also notice that the correlator C(k±) obeys the sum rule 

^C(k x ) = e-*<°> = 1 - (83) 

Numerical results as well as the asymptotic formula for C(kj_) are displayed 
in figure 0. 
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Figure 4: Value of C(k±) as a function of k±/A QCD for Q S /A QCD = 10. The 
open circles are the result of a numerical evaluation, and the solid line is the 
asymptotic formula given in Eq. (p2|). The shaded area indicates the region 
where k± < A QCD . One can see on the numerical curve the onset of saturation 
at low values of kj_. In the upper right corner, we have plotted the same data for 
the quantity k 4 L C(k±)/2Ql, which allows for an easier reading for large values 
of kj_. 
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